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Abstract 

The behaviour of the quark coefficient functions in deep-inelastic scattering is investigated 
for large values of the Bjorken variable x. By combining results of soft-gluon resummation 
and fixed-order calculations, we determine the coefficients of the four leading large-x 
logarithms, a k s [{hi(l — x)} 2k ~ l /(l — x)] + , I — 1, . . .4, to all orders in the strong coupling 
constant a s . This result includes two more terms for the three- loop coefficient functions 
than previously specified in the literature. The effect of the fifth logarithmic contribution 
is approximately evaluated. The terms derived here are required, but also seem to be 
sufficient, for a reliable representation of the coefficient functions at large x. 



Knowledge of the proton's quark distributions at large momentum fractions x is relevant 
to the search for new phenomena at high-energy colliders like HERA, Tevatron and 
the future LHC These distributions are usually inferred from data on structure functions 
in deep-inelastic scattering (DIS). For instance, data even up to x = 0.98 have recently 
been employed |1J to scrutinize a possible non-standard behaviour of the large- a; quark 
distributions 0] suggested in connection with the initial HERA high-Q 2 anomaly 0. In 
perturbative QCD the link between the quark distributions and the structure functions 
Fi(x, Q 2 ), i — 1, 2, 3, is given by the coefficient functions^] 

C i>q (x,Q 2 ) = 5(l-x) + J2a k s ^(x) (1) 
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with a s = a s (Q 2 )/(4ir). The expansion coefficients c\J(x) and their Mellin-iV space coun- 
terparts c\J contain universal (i.e., i-independent) soft-gluon contributions of the form 

(In' N + subleading terms) , (2) 
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1 = 1,..., 2k. At sufficiently large x / large N, these terms spoil the convergence of finite- 
order approximations to Cj i9 and need to be resummed to all orders. It is thus important 
to derive a reliable estimate of their impact at the x- values included in data analyses. 

The soft-gluon resummation has been carried out for the leading and next-to-leading 
logarithms (in the sense of eq. (|8D below) some time ago @, ||. In this letter we combine 
these results with terms from the one- and two-loop coefficient functions c\ ' (x) calcu- 
lated in refs. || [/J. This information fixes the first four towers, j = 0, 1, 2 and 3, of large-x 
logarithms, i.e., the coefficients of In 2 ~ 3 N in c\J at each order k. It also facilitates an 
estimate of the effect of the fifth tower, j = 4. As shown below, these four to five towers 
are required, but also seem to be sufficient, for a realistic estimate of the higher-order 
effects at large x. 

Up to terms which vanish for N —* oo, the iV-space coefficient functions take the form 

C^Q 2 ) = (1 + a s g 01 + a% 2 + ...)• exp [G N (Q 2 )} . (3) 

The first factor collects contributions which are constant for N — * oo. Adopting the MS 
scheme, the first-order term reads || 

g 01 = (-9-2C 2 + 2 7 2 + 3 7e )CV • (4) 

Here d stands for Riemann's ^-function, 7e is the Euler-Mascheroni constant, and C F = 
4/3 in QCD. The corresponding second-order term can be readily inferred from ref. 0. 

lr rhe renormalization and factorization scales are chosen as \j? t = [il = Q 2 throughout this paper. The 
additional terms arising for other choices can be deduced from renormalization- group constraints. 



The function G N (Q 2 ) is the object of the soft-gluon resummation. It is given by 0, 
G N (Q 2 ) - 
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+ 0(a s (a s \nN) k ) 
with A(a s ) = 4 C F a s + 8 C F K a 2 + . . ., B = -6 C F a s + . . . and 
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Ca = 3 in QCD, and Nf denotes the number of effectively massless quark flavours. These 
terms of A(a s ) and B(a s ), together with the lowest two coefficients of the /3-function, 



-C A - 2 -N f , 
3 A 3 f ' 



are sufficient for deriving the next-to-leading logarithms. 
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The integrals occurring in eq. (|5]) can be found in ref. ||. Using the abbreviations 
L = IniV and A = f3oa s L, the result can be written as 

G N {Q 2 ) = L gi {a s L) + g 2 {a s L) + a s g 3 {a s L) + ... (8) 



with 
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(9) 



A + ln(l-A) 



A + ln(l-A) + -ln 2 (l-A) 
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Here O^j = 1 for k > j and #£,,■ = else, and Si(k) = J2j=i 1/j- The function g 3 in eq. (§), 

<? 3 (a s L) = ^^(a^)*" 1 . ( n ) 
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is presently unknown except for its leading term which can be determined from ref. 
Its coefficient g 32 reads 

/3155 , n . 22 A oA 22 2 367 \ _ „ 
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We are now ready to evaluate the resummed large- N coefficient function (|3|). Expan- 
sion of the exponential yields 



CUQ 2 ) = ! + £«} (c kl L* k + c k2 L^ + c k3 L™-> + c M L^ + 0(L 
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Inspection of eq. (|14|) reveals that the coefficients up to c fc4 are indeed fixed by eqs. (f|), 
(0), ([U]) and (|T2|).f| Especially the four leading large-x terms of the three-loop (k = 3) 



coefficient functions q are thus determined. Transformation to x-space || yields 
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The first two terms of eq. (|15 ) have already been given in ref. [11]]. Besides as a useful 



.(3), 



cross check for a future exact calculation of c\'(x), the above result can also be used for 



improved approximate reconstructions |IIJ of this function along the lines of refs. [D], [12|] 



2 The agreement of the predictions rtlj) for c 21 , c 22 and c 23 with the results of the full 0{a 2 s ) calcula- 
tion M provides a non-trivial check of the resummation (P)-([lO|) of the next-to-leading logarithms. 



The corresponding higher-order results are straightforward if cumbersome. The numerical 
values of the coefficients c kl in eq. fllil) are presented in table 1 for Nf = 4 and k < 10. 
In contrast to g 2 in eq. ( |TU| ) which exhibits a pole at N = exp[l/(/3 a s )], the sum flT3] ) 
converges for all A" as long as the number of towers included is finite. 
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Cfcl 


C fc2 


C k3 


C k4 


C k5 


1 


2.66667 


7.0785 











2 


3.55556 


26.2834 


40.760 


-67.13 





3 


3.16049 


44.9210 


238.885 


470.82 


#33 - 1235.23 


4 


2.10700 


48.7090 


477.854 


2429.46 


cn033 + 3600.12 


5 


1.12373 


38.3254 


581.518 


5015.18 


On #33 + 22963.9 


6 


0.49944 


23.5617 


505.972 


6432.95 


C31 9 33 + 50185.0 


7 


0.19026 


11.8592 


340.954 


5933.61 


C41 g 33 + 67307.5 


8 


0.06342 


5.0464 


186.822 


4244.45 


C51 933 + 64858.9 


9 


0.01879 


1.8583 


86.041 


2467.72 


Oil #33 + 48498.6 


10 


0.00501 


0.6028 


34.118 


1204.34 


en #33 + 29487.7 



Table 1: Numerical values of the coefficients c kl in eq. (0) for A 7 / = 4. Also shown are 
the corresponding (incomplete) results c k5 for the fifth tower of logarithms in eq. (|i~3|) . 



Before turning to the effect of these higher-order contributions, it is instructive to 
compare the a 2 s and oP s parts of eq. (|13|) with the respective exact and approximate 
results for cjji and c 2 i [11], P^1-P| For this purpose the coefficient function is convoluted 
(for as = 0.2, a value typical for scales probed in fixed-target DIS) with a simple, but 
characteristic model f(x) of the large- a: quark distributions. The second-order comparison 
is shown in fig. 1, its third-order counterpart in fig. 2. Two and three terms in the large- A" 

expansion flT3"|) (i.e., In Amandin A 7 at two-loop and in N, In N, and In A^ at three-loop) 

(3) (3) 

turn out to be sufficient for good approximations to q' <8> / and c\' (g> /, respectively, 

at large x. It should be noted that it is essential for this fast convergence to use the 

(2) 
expansion in A^-space. For instance, c 2 ' ® / is severely overestimated if only the two 

leading x-space terms, [ln 3 (l — x)/(l — x)] + and [ln 2 (l — x)/(l — x)]+, are kept [ll|. In 



the all-order case, the problem of the x-space expansion has been elucidated in ref. fI3 



The higher-order soft-gluon corrections, k > 4 in eq. (|T3p , are illustrated in the same 
manner in fig. 3. Under these conditions the terms up to k = 8 are sufficient for an accurate 
representation up to z = 0.99. As above the Mellin inversions are performed using the 

3 These approximations are based on the first five even- integer moments of c 2 i(x) UM, supplemented 
by the first two coefficients in eq. (|l5|). c 2 i(x) is rather tightly constrained at x > 0.5 by this information. 



standard contour |L5] (see also the discussion of the infinite-order limit in ref. [ |i"4|| ) 

I f PC rc+(i-S)oo 1 

A W = T~- / +/ diVx-^Ar 
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— I dz Im 
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exp(z0) x- c " zexp ^) ^ =c+2exp(w J (16) 



using c = 2 and <p — 3/4 7r corresponding to 5 = 1. The fourth tower of logarithms, 
J2k c k4 a> s l 1 * 2 ~ 3 ^, yields a very large contribution, e.g., it exceeds the effect of the first 
three towers by a factor of 1.5 at x = 0.9. In order to assess the convergence of the 
large- N expansion, it is useful to notice that the fifth tower, J2k c k5 a s hi 2fc_4 iV, is fixed 
by available information (including the parameter g Q2 in eq. (f|)), except for the second 
term g 33 (a s ln N) 2 of the function g 3 in eq. (§). As shown in table 1, the impact of the 
unknown constant g 33 (of which the order of magnitude can be inferred from h — 3) is 
small at k > 5. The effect of the fifth tower can thus be estimated, and turns out to 
be rather moderate, e.g., about 25% at x = 0.9. This indicates that four to five towers 
are sufficient for a valid estimate of the soft-gluon corrections to the quark coefficient 
functions at large x. 

To summarize, we have determined the four leading towers of large-a; / large- iV soft- 
gluon logarithms entering the quark coefficient functions in DIS, and estimated the fifth 
tower. There is evidence that - provided the resummation is applied in iV-space, avoiding 
uncontrolled lower In TV terms due to the Mellin transform (0) - these four to five towers 
are sufficient for a reliable estimate of the higher-order contributions at large x. These 
terms should thus be included in future analyses of structure function data at very large x. 
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Figure 1: Convolution of the second-order contribution to the coefficient function C^ q 
with a typical input shape. Shown are the full result [|7| for c 2 l and four large- N approx- 
imations, in which the terms ~ In N, In N, In N and IniV in eq. ([T^) are successively 
included (e.g., the long-dashed curve includes the ln 4 iV and In 3 N parts). 
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Figure 2: Convolution of the third-order term c^' q of the coefficient function C2, q with 
a typical input shape. The four leading large-iV approximations are compared with a 
parametrisation [DI of the full result based on the lowest five even-integer moments [II 



The spread of the full curves indicates the residual uncertainty of this parametrisation. 
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Figure 3: The effect of the higher-order soft-gluon corrections, Y,k=4 Sj=o c kj a s ^ ~ 3 N \ 
to the quark coefficient functions at large N. The results including up to four towers 
(jm — 3) are derived from the exact coefficients in eq. (0). The impact of the fifth tower 
(jm = 4) is estimated using g 33 



500 for the coefficients c k5 at k > 5 in table 1. 



